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■ Abstract 

x 

' An analysis of electromagnetic corrections to the (dominant) octet 

K — > 7T7T hamiltonian using chiral perturbation theory is carried out. 
Relative shifts in amplitudes at the several per cent level are found. 



1 Introduction 



In this paper, we present a formal analysis of electromagnetic (EM) radiative 
corrections to K — ► tttt transitions.^ Only EM corrections to the dominant 
octet nonleptonic hamiltonian are considered. Such corrections modify not 
only the original AI = 1/2 amplitude but also induce AI = 3/2,5/2 con- 
tributions as well. By the standards of particle physics, this subject is very 
old Yet, there exists in the literature no satisfactory theoretical treat- 
ment. This is due largely to complications of the strong interactions at low 
energy. Fortunately, the modern machinary of the Standard Model, espe- 
cially the method of chiral lagrangians, provides the means to perform an 
analysis which is both correct and structurally complete. That doing so 
requires no fewer than eight distinct chiral langrangians is an indication of 
the complexity of the undertaking. 

There is, however, a problem with the usual chiral lagrangian methodol- 
ogy. The cost of implementing its calculational scheme is the introduction of 
many unknown constants, the finite counterterms associated with the regu- 
larization of divergent contributions. As regards EM corrections to nonlep- 
tonic kaon decay, it is impractical to presume that these many unknowns 
will be inferred phenomenologically in the reasonably near future, or perhaps 
ever. As a consequence, in order to obtain an acceptable phenomenologi- 
cal description, it will be necessary to proceed beyond the confines of strict 
chiral perturbation theory. In a previous publication ||], we succeeded in 
accomplishing this task in a limited context, K + — > tt + -k° decay in the chiral 
limit. We shall extend this work to a full phenomenological treatment of 
the K — ► 7T7T decays in the next paper [Q] of this series. 

The proper formal analysis, which is the subject of this paper, begins 
in Sect. 2 where we briefly describe the construction of K — > tttt decay 
amplitudes in the presence of electromagnetic corrections. In Section 3, 
we begin to implement the chiral program by specifying the collection of 
strong and electroweak chiral lagrangians which bear on our analysis. The 
calculation of K — * tttt decay amplitudes is covered in Section 4 and our 
concluding remarks appear in Section 5. 



x We restrict our attention to EM corrections only and omit consideration of m u ^ m,d- 
See however Ref. [l| 



1 



^< -^x 

(a) (b) (c) 

Figure 1: Some electromagnetic contributions. 

2 Electromagnetism and the K — > 7T7t Amplitudes 

There are three physical .fT — > 7T7r decay amplitudes ,0 

^lft-0^ 7r + 7r - = , ^4x ^7r°7r — -4-00 , -^ftT+^vr+Tr = -4+0 • (1) 

We consider first these amplitudes in the limit of exact isospin symmetry 
and then identify which modifications must occur in the presence of electro- 
magnetism. 

In the 7 = 0,2 two-pion isospin basis, it follows from the unitarity con- 
straint that 

A+- = A e iS ° + ^A 2 e iS * , 
A 00 = A e iSo - V2A 2 e iS2 , (2) 
A +0 = \A 2 e^ . 

The phases 5q and 5 2 are just the 1 = 0,2 pion-pion scattering phase shifts 
(Watson's theorem), and in a CP-invariant world the moduli Aq and A 2 are 
real- valued. The large ratio Aq/A 2 ~ 22 is associated with the AI =1/2 
rule. 

When electromagnetism is turned on, several new features appear: 

1. Charged external legs experience mass shifts (c/Fig. |l|(a)). 

2. Photon emission (c/Fig. |l](b) ) occurs off charged external legs. This 
effect is crucial to the cancelation of infrared singularities. 

3. Final state coulomb rescattering (c/Fig. |l|(c)) occurs in K° — > ir + ir~ . 
2 The invariant amplitude A is defined via out(7r7r|-K}in = i(2ty) 4 S^ (p Q ut — Pin) (iA). 



2 



7 




Figure 2: Leading electromagnetic correction to K — > 7T7t. 



4. There are structure-dependent hadronic effects, hidden in Fig. 1 within 
the large dark vertices. In this paper, we consider the leading contri- 
butions (see Fig. |||) which arise from corrections to the AI = 1/2 
hamiltonian. 

5. There will be modifications of the isospin symmetric unitarity relations 
and thus extensions of Watson's theorem. 

Any successful explanation of EM corrections to K — * tttt decays must ac- 
count for all these items. 

An analysis || of the unitarity constraint which allows for the presence 
of electromagnetism yields 

A + _ = (A + 5A^ m ) e^o+To) + _L (a 2 + $Af) e { ^ + ^ , 

v2 

„4 00 = (A + 5Al m )e^ 6o+ ^ -V2(A 2 + 5Al m )e i ^ + ^ , (3) 



Ao = l(A 2 + 6Af 



em\ J{$2+l' 2 ) 



to be compared with the isospin invariant expressions in Eq. (|2|). This 
parameterization holds for the IR-finite amplitudes, whose proper definition 
is discussed later in Sect. 4.3. Observe that the shifts <L4^~ em and 72 in 

.4+0 are distinct from the corresponding shifts in A-\ and Aqq. This is a 

consequence of a AI = 5/2 component induced by electromagnetism. In 
particular, the AI = 5/2 signal can be recovered via 



A 

A/2 - — 



A+- - Aoo - V2A + o 



(4) 
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3 Chiral Lagrangians 



The preceding section has dealt with aspects of the K — > tttt decays which 
are free of hadronic complexities. In this section and the next, we use 
chiral methods to address these structure-dependent contributions. The 
implementation of chiral symmetry via the use of chiral lagrangians provides 
a logically consistent framework for carrying out a perturbative analysis. 

In chiral perturbation theory, the perturbative quantities of smallness 
are the momentum scale p 2 and the mass scale x = ^Bqui, where m is the 
quark mass matrix. In addition, we work to first order in the electromagnetic 
fine structure constant a, 

A = if ) +aif ) + ... . (5) 

Our goal is to determine the 0(a) components aA.^ ■ The fine structure 
constant thus represents a second perturbative parameter, and we consider 
contributions of chiral orders O(e 2 p ) and 0(e 2 p 2 ), 

aA? = A (*V) +A (eV) . ( 6 ) 

We shall restrict our attention to just the leading electromagnetic cor- 
rections to the K — > tttt amplitudes. Since the weak AI = 1/2 amplitude is 
very much larger than the AI = 3/2 amplitude, our approach is to consider 
only electromagnetic corrections to AI = 1/2 amplitudes. As a class these 
arise via processes contained in Fig. ||, where gs is the octet weak coupling 
defined below in Eq. (|l3|). 

We adopt standard usage in our chiral analysis, taking the matrix U of 
light pseudoscalar fields and its covariant derivative D^U as 

U = expia^k/Fn) {k = 1, . . . , 8) , DfJJ = d^U + ie[Q, U]A^ , (7) 

where Q = diag (2/3,-1/3,-1/3) is the quark charge matrix and is the 
photon field. The remainder of this section summarizes the eight distinct 
effective lagrangians (strong, electromagnetic, weak and electroweak) needed 
in the analysis. 

3.1 Strong and Electromagnetic Lagrangians 

In the AS" = sector, we shall employ the strong/electromagnetic lagrangian 

42 = ^ Tr (D,UD»W) + 3l Tr ( x tft + ^t) , ( 8) 
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(2) 

where F is the pseudoscalar meson decay constant in lowest order. will 
be used to produce O(e p 2 ) and 0{e 1 p 1 ) vertices in our calculation. 

(2) 

The lagrangian will generate (via tadpole diagrams) strong self- 
energy effects on the external legs in the K — » irir transitions. In order 
to regularize these divergent contributions, one employs the lagrangian || 
£gtr- It is not necessary to write out this well-known set of operators, but 
simply to point out that the resulting wave function renormalization factors 
Z n and Zk obey 

1 Z^^fZrK 



F 2 F K F* ' (9) 

up to logarithms. This explains the presence of F 2 Fk in formulae such as 
Eqs. (H),© in Section 4. 

Two other nonweak effective lagrangians enter the calculation. The first 
is associated with electromagnetic effects at chiral order O(e 2 p ), 

4% = g eim Tr (QUQrf) , (10) 

where the coupling g ems is fixed (in lowest chiral order) from the pion elec- 
tromagnetic mass splitting, 

9ems = ^-SM 2 . (11) 

The second extends the description to chiral order 0(e 2 p 2 ). We need only 
the following subset of the lagrangian given in Ref. p], 



C (2) = p2 2 



K\ 



Tr (DMD^) ■ TiQ 



\2 



+ K2 Tr (D^UD^) ■ Tr [QUQU^ 
+ k 3 (Tr (D^QU) ■ Tr (d^QU) 

+ Tr (D^UQU^ • Tr (d»UQU^ 
+ K 4 Tr (DprfQU) ■ Tr (D^UQU ] ^j (12) 
+ k 5 (Tr (p^DfJJQ) + Tr (d^UD^Q)) 

+ kq Tr (D^D^UQU^QU + D^UD^QUQU^ 



Although the finite parts of the coefficients k±, . . . ,Kq remain unconstrained, 
see however Refs. [0, ||, 10 [ for model determinations. 
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3.2 Weak Lagrangians 

The \ AS\ = 1 octet lagrangian begins at chiral order p 2 , 

4 2) = 58 Tr (\ 6 D„UD»Ut) , (13) 



with gg ~ 6.7 • 10~ 8 F 2 fit |jl6| from K — > -kit decay rates. We use this to 
generate O(e p 2 ), 0(e 1 p 1 ) and O(e 2 p ) vertices. 

Two chiral lagrangians will serve to provide counterterms for removing 
divergent contributions. The first [11] is the octet \AS\ = 1 lagrangian at 
chiral order p 4 , 

4 4) = N 5 TrX 6 (u X j + xU ] ) dJJ&rf + d^Ud^tf (u X ] + xU ] ) 
+ N 6 Tr X e Ud^ ■ Tr (x}d»U - x^C/ f ) 
+ N 7 Tr A 6 (U X ] + X^ f ) • Tr d^Ud^tf 
+ N 8 Tr Xad^Ud^tf ■ Tr (u ] X + X f ^) 

+ N 9 Tr A 6 d^Ud^U^ (x^ 1 " - U X ] ) - [xU ] - U X ] ) d^Ud^U ] 
+ N w Tr A 6 W + X UhU ] + U^X^) 
+ N n Tr A 6 (yx 1 + X rf) • Tr (rf x + X ] u) 
+ N 12 Tr A 6 (c/ X W + xU i xU 1 _ c/^X^ 1 ) 
+ AT 13 Tr A 6 (c/ X f - Xf/ f ) • Tr (c/ X f - X^) . (14) 
At present, little is known of the finite parts of the couplings {Nk}. 

3.3 Electroweak Lagrangians 

The \AS\ = 1 lagrangian at chiral order O(e 2 p ) is 

4Sw=5emwTr (WQETf) , (15) 

where g emw is an o priori unknown coupling constant. It has been calculated 
recently in Ref. ||, 

W = (-0.62 ±0.19) g 8 5M 2 . (16) 
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We note in passing that despite the presence of just one charge matrix Q 
the lagrangian of Eq. (|D|) indeed describes 0(e 2 ) effects. A second factor of 
Q could be decomposed into a combination of the unit matrix and the 3x3 
matrix Q = diag (1,0,0). The contribution from Q would vanish, leaving 
the form of Eq. flj~5[ ) . 

The second operator that we use to provide counterterm contributions 
is the |AS| = 1 lagrangian at chiral order 0(e 2 p 2 ). In terms of the notation 
Lu = illdnW, we have 



=e 2- 

""emu) ^ 



si TrA 6 [Q,V2L^] + 
+,s 2 Tr A 6 (QUQrtLpW + L^UQU^Q) 

+s 3 Tr A 6 [Q, L f ,UQU^L^]+ + s A Tr A 6 UQU^}+ ■ Tr UQlfllf 
+s 5 Tr A 6 (QUQU^ X U ] + Ux ] UQWQ) 

+s 6 Tv\ G [x,U^} + • Tr UQWQ (17) 
+s 7 Tr A 6 

+s 8 Tr (\^Ud^) ■ TrQ 2 

+s 9 Tr (Xed^Ud^U^ • Tr UQU^Q . 



The first six operators in the above list appear in Ref. |12|| . The remaining 
three are also required for our analysis. To our knowledge, none of the 
divergent or finite parts of the {s n } are yet known. 



4 Calculation of Leading EM Corrections 

The leading EM corrections arise from the processes of Fig. [l] and Fig. ||. 
Contributions to Fig. |2| occur in two distinct classes, those explicitly con- 
taining virtual photons (Fig. |3|) and those with no explicit virtual photons 
(Fig. [|). The latter are induced by EM mass corrections and by insertions 
of 

<?emw I n Figs. |3[^|, the larger bold- face vertices are where the weak inter- 
action occurs. 

The integrals which occur in our chiral analysis are standard and already 
appear in the literature (e.g. see Ref. H or Ref. 0). It suffices here to 
point out that all divergent parts of the one-loop integrals are ultimately 
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(a) (6) (c) (d) 

Figure 3: Explicit photon contributions in K + — > 7r + 7r°. 



expressible in terms of the ci-dimensional integral 
1 



A(M 



dk 



d-4 



k 2 -M 2 



2 - iM 2 , M 2 \ 
^A-^log l-f U... 



(18) 

where dk = d k/(27t) d is the integration measure, [i is the scale associated 
with dimensional regularization and A is the singular quantity 



A = 



1 



16vr 2 



^-"(log^-T + l) 



(19) 



Each amplitude in the discussion to follow will be expressed as a sum of a 
finite contribution and a singular term containing A. 



4.1 Summary of 0(e 2 ) Amplitudes 
We begin with the O(e 2 p ) amplitudes, 

V2 



A 



(e 2 p°) 



FkF 2 



(gsSM 2 



+ 9e 



' -^00 



, A 



( e 2 p 0) _ A + 

+0 



(eV) 



V2 



(20) 

Although these have already been determined in Ref. ||, we include them 
here for the sake of completeness. They are finite-valued and require no 
regularization procedure. 

Next come the amplitudes of order e 2 p 2 , expressed as 



^(eV) _ ^(cxpl) _j_ ^(impl) + ^(ct) 



(21) 



The superscript 'expl' refers to Figs. [11(a), (c) and Fig. [| where virtual pho- 
tons are explicitly present, whereas superscript 'impl' refers to Fig. |3| where 
EM effects are implicitly present via EM mass splittings and g e mw insertions. 
The final term A^ is the counterterm amplitude. 
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4.1.1 Diagrams with Explicit Photons 

We turn first to the class A^ expl ^ of explicit photonic diagrams. For these 
contributions, it is consistent to take meson masses in the isospin limit. We 
find 

FkFI 



V2g 8 



^(expl) = ^ M 2 K _ M 2^ . aB+ {mi) 



a 

AlT 



7Mf 



3M£ 




6// 



F KF n .(expl) 

V2 9s A " 
FkFI 



(22) 



A 



(expl) 
+0 



a 

Att' 



-Mi 



7-3 hi- 



Ml 



Qn d ~ i e 2 Ml\ 



The quantity , which appears in the above expression for .A^* pl \ 

is associated with the processes of Figs. |](a),(c). Due to such processes, 
the weak decay amplitudes A% will develop infrared (IR) singularities in the 
presence of electromagnetism. To tame such behavior, an IR regulator is 
introduced and appears as a parameter in the amplitudes. For our work, 
this takes the form of a photon squared-mass m 2 . is given by 



B + -{m: 



1 

Air 



M 

2o(/3) In -f 



m 



2 1 + 2 1 + 



7 



26 



1-/3 



+ in 



1 + 2 





In 



M 2 K 2 



7714 







where 



and 



a(0) 
HP) 

m 



i + 







i + ^ ln i 



(1- AM 2 /Ml-) 1/2 
-0 



(23) 
(24) 



26 



1 + 



7T + m in ■ 



1-0 A0- 



2 



1 + 

20 



2/ 



0-1 

20 



(25) 



dt - In 1 1 

t 1 



t 



Notice that the function B^ is complex, and both its real and imaginary 

parts have a logarithmic singularity as 777 7 —* 0. The solution to this problem 
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(a) (6) (c) (d) 

Figure 4: Diagrams without explicit photon contributions. 



is well known; in order to get an infrared-finite decay rate, one has to consider 
the process with emission of soft real photons, whose singularity will cancel 
the one coming from soft virtual photons. We shall be more explicit on this 
point in Sect. 4.3. 

The amplitudes A+^F^ and ^4,^q P ^ each contain an additive divergent 
term (proportional to A) and also depend on the arbitrary scale jx introduced 
in dimensional regularization of loop integrals. Both these features will 
require the introduction of counterterms. 



4.1.2 Diagrams without Exploit Photons 

Next comes the class _4( im P [ ) of dia grams in Fig. || not containing explicit 
photons. For such contributions, one must be sure to include all possible 
effects of chiral order O(e 2 p ) and 0(e 2 p 2 ) and treat the various terms in 
a consistent manner. Thus for the contributions to Fig. ^j, isospin-invariant 
meson masses are used in amplitudes involving x C^j. and x 

(2) 

C 8 , whereas electromagnetic mass splittings appear in amplitudes involving 

(2) (2) 

£g tr x £g . We write the results as sums of complex- valued finite amplitudes 
Fiifi) and divergent parts, essentially the amplitudes T>i, 

A {imp\) = Ue r ^ + r ^ + ^d-4 p . ( • = + _ 0Q) +Q) _ (2g) 

The scale-dependence in J-i(^) comes entirely from its real part IZe Tii/J.). 
We express the Tie Ti in terms of dimensionless amplitudes a\ impl \ 

Ke W = Vi^af^M , (27) 

with r/_| = 7]qq = \/2, V+o = 1- Since the a[ irapl \fi) coefficients have rather 

cumbersome analytic forms, we find it most convenient to express them in 
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the compact form 

4 mp V) 

where 



(M)5M n (g) g 



9 



(M)jAg (g) g 



In 



1 GeV 



9emw/98 ■ 

The coefficients appearing in Eq. (p8|) are given in Table 1. 

Table 1: Values of Coefficients in Eq. (|29|) 



(28) 
(29) 





6 (M) 


6 W 


(M) 


Jg) 


i = +- 


0.0160 


-0.0409 


-0.0078 


-0.0445 


i = 00 


-0.0170 


-0.0224 


-0.0371 


-0.0176 


i = +0 


-0.0265 


-0.0220 


-0.0419 


-0.0357 



The finite functions also have imaginary parts Im T% which arise entirely 
from the processes in Fig. ||(c). From direct calculation we find 



F K F^ 
V2g 8 



Im J-+- 



JL 

16vr 



Ml 



5M 2 +g 



-\^-2\{Ml-Ml)ml 



ft 2 



2 ip2 



F K FjF 
V2g 8 



■Im T{ 



oo 



2 772 



F K FjF 
98 



■Im J-. 



+o 



JL 

16vr 



JL 

32vr 



Ml 



5Ml + g 



+ 2- 



M 2 K - Ml 5Ml 



P 2 



(30) 



M 



K 



2M 2 



(5M 2 



9 



where j3 is defined in Eq. (24). As a check on our calculation, we have 
verified that the above results are identical to those obtained from unitarity. 



The singular parts of A 
F 2 F K Fl 



(impl) 



are embodied by the D-functions, 



V2g 8 
F 2 F K Fl 

V2g 8 
F 2 F K Fl 

98 



V+- 
^oo -- 



M 2 K 



2 n 2 



+ Mi 



(M 



Ml 



™dM 2 + 3g 



M 



K 



—8 Ml H g 

3 71 18 



+ Mi 



105M 2 + 7g 



,.„, 86 
A8M 2 + -g 



(31) 
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(a) (b) ~ (c) 

Figure 5: Counterterm contributions. 

To arrive at the above, we have used both the correspondence between 5M 2 
and <7 ems given in Eq. (|TTj) and also the relation 



M 2 ± -M^ = M 2 K+ -M 2 R0 , (32) 
in the evaluation of loop integrals. The latter follows from Dashen's theo- 



rem [15] and is justified since terms violating Dashen's theorem would begin 
to contribute at the higher chiral order e 2 p 4 . 

4.2 The Regularization Procedure 

In order to cancel the singular A-dependence in the K — > tttt amplitudes, 
it is necessary to calculate all possible counterterm amplitudes which can 
contribute. These enter in a variety of ways, as shown in Fig. || where the 
small bold-face square denotes the counterterm vertex. For Figs. |5|(a),(b) 
the counterterm vertex has \AS\ = 1 whereas in Fig. ||](c) it has AS = 0. 

4.2.1 Counterterm Amplitudes 

Using the lagrangians ci^Lv and £ems we determine the counterterm 
amplitudes to be 

F 2 F K F 2 & = 
V2g s + 

2 /'jp2/V. ATT 8 TT\ , XH4-2I 



m k WF\Xx - AUi - -U 2 ) + 5M 2 (8N 7 - 4N 8 - 4iV 9 ) 
+ Ml (e 2 F 2 (X 2 + AUt + -U 2 ) - SM 2 (4N 5 + 8N 7 + 2N 8 ) 



F 2 F K F 2 A {ct) _ ( nf2 _ n/r 2\ 2t?2 v.. _ \n _ : _ -nr. (33) 



^00 - m - \u 2 - 2U 3 



F2FkF ' A^ < > = Ml ( e 2 F*X 3 - SMl(iN ; , + 4JV 8 
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+ Ml ( e 2 F 2 X 4 - 5M'l(2N 8 + 47V 9 ) 



where the {N} are coefficients in the \AS\ = 1 lagrangian 

4 4) of Eq. fllj, 

(2) 

the {Ui} are combinations of coefficients in the AS = lag rangian ^Cems of 
Eq. (0), 

U\ = Ki + K 2 j #2 = ^5 + K 6 i ^3 = -2^3 + «4 , (34) 

and the {Aj} are combinations of coefficients in the |A5"| = 1 lagrangian 
C%L of Eq. (pi), 



4 1 2 2 ^ 2 

Xl = ~-S! - -S 2 + -S 3 + -S 5 - 4s e + ~S 7 + Sg + Sg , 

4 2 4 4 2 
X 2 = ~si - -s 2 + -s 3 + -s 5 + 4s 6 - -s 7 - s 8 - s g , 
y y y o o 

2 14 2 2 
^3 = -gSi - g«2 + g«4 + gS5 + gS 7 , (35) 

2 2 4 4 2 
X 4 = - ai + - S8 - r 4+ 3.5-3.7 , 

2 2 

-^00 = g («1 + «2 + S3) + 3S4 + S 8 + Sg , 

4.2.2 Removal of Divergences 

The counterterms themselves have finite and singular parts, 



(36) 



The coefficients rii,Ui of the divergent parts of Ni,Ui have already been 
specified in the literature JO], [?J and hence the //-dependences of N- r \ 
are known from the renormalization group equations. We infer the Xj co- 
efficients in this paper by canceling divergences in the 0(e 2 p 2 ) amplitudes. 
Upon combining results obtained thus far, we find the new results 

ISM 2 g 

X 00 — —7; 9 no — >J 



Ni 


= riiH d - 




Ui 




~ A \ + ul r) 


Xi 


= xni d - 


- 4 x + xl r) 



x 1 = 3 + 



3e 2 F 2 " e 2 F 2 ' 

27 5M 2 13 g 



2 e 2 F 2 2 e 2 F 2 
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X-2 



X ?> 



3-18 



5Ml 



e 2 F 2 

7 6Ml 89 



9 

e 2 F 2 

9 



(37) 



3 e 2 F 2 18 e 2 F 2 
n 5M 2 86 g 



e 2 F 2 g e 2 F 2 



where we recall g = g e mw/#8- 



4.3 Removal of Infrared Singularities 

Removal of the infrared divergence from the expression for the decay rate is 
achieved by taking into account the process K° — ► 7r + 7r~(n7). For soft pho- 
tons, whose energy is below the detector resolution u, this process cannot be 
experimentally distinguished from K° — » ir + ir~, so the observable quantity 
involves the inclusive sum over the h + it~ and 7r + 7r~(n7) final states. 

At the order we are working, it is sufficient to consider just the emission 
of a single photon. The amplitude for the radiative decay is given in lowest 
order by 



A 



+-7 



FkF 2K ' 



Mi 



e-p + e-p. 



q-p + q-p. 



(38) 



where e and q are the polarization and momentum of the emitted photon. 
The infrared-finite observable decay rate is 



where 



r+- 7 (w) 



2M K 
1 

2M K Je. 



,<UJ 



7 I *^^H — 7 1 



(39) 

(40) 
(41) 



and is the differential phase space factor for each process. The infrared 
divergent (IRD) part of r_| is seen to be 



,(IRD) 



1 



2M K 



FkF 2 



(M 2 K -Ml 



2a^e5 + _(m 7 ) . (42) 



Equation (42) displays explicitly the singularity and shows that the imag- 
inary part of -B-i (jrij) has no observable effect at this order. This result 
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has been shown to be true to all orders in a Jl7| , S3]. For T H 
the following expression, up to terms of order lo/Mk, 



(lo) we get 



2M K 



V2g 8 (m2 



FkFI 



K 



Mi 



I-\ (m 7 ,u;) 



where 



with 



/-l (m 7 , w) 



Q: 
7T 



(43) 



(44) 



1 1 + 1 + /3 



2 r 



/5 ln i 



+ 



2/3 



2/(-/3) - 2/09) + / 



1 + /? 



/3 



+ I In ln(l - /3 2 ) + ln2 In - 

2 1-/3 v ^ ' 1 + /3 



.(45) 



From these explicit expressions of B^ (lo) and I-\ (m 7 ,u) it is easy to 

see that the combination 2alZeB^ (m^) + 1-\ (m 7 ,u;) does not depend on 

the infrared regulator m 7 . However, this combination has a dependence on 
the experimental resolution lo. To obtain a meaningful prediction therefore 
requires knowledge of the experimental treatment of soft photons. A careful 
discussion of this point will appear in Ref. ||. 

A generalization of the above considerations beyond the order 0(e 2 p 2 ) 
in ChPT leads to the following parameterization, 



i 



[to 



A ( ?l + a A 



(1) |2 



2M K 

where to first order in a, 

G-i (lo) = 1 + 2alZeBj l (m 7 ) + I-\ (m 7 ,a;) 



(46) 



(47) 



With the prescription of dropping the term proportional to B^ in the 

photonic loop contribution, the electromagnetic amplitude aA^P can be 



read from Eqs. <m,&M),(M)- 



4.4 The Finite Amplitudes 

The physical amplitudes will be complex- valued functions, as dictated by 
unitarity. The real parts are obtained by combining the finite loop ampli- 



tudes (Eq. Q22) for ^J cxpl) and Eqs. (E7]),(| 



along with Table 1 for A^ ) 
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with the counterterm amplitudes of Eq. (|3 



Tie A, 



(eV) 



g&M j 



K r 



F%F K 



Tie af oop) + 



In order to make the scale-dependence of IZe af oop ^ explicit, we write 



Tie a) 



(loop) 



bi + Ci In 



1 GeV 



(48) 



(49) 



Numerical determination of the above quantities will depend on g$ (obtained 
from Ref. [jlgl), 5M% and ij emw (given in Eq. (|l6|) ). We obtain the central 
values 



b + - 
boo ■- 
b+o 



11.8 • 10~ 3 
-0.5 • lO- 3 . 
-1.3 • 10" 3 



coo z 
c+o 



7.1 • 10" 3 , 
-3.9 • 10-3 , 
-2.7 • 10 -3 



(50) 



The imaginary parts of the physical amplitudes can be either determined 
from unitarity or read off from Eqs. (|26|) , (|3(i|) . Of most interest is the EM 
shift in A2, as only it receives the A0/A2 (A/ = 1/2) enhancement, 



S(Im A e 2 m ) 



JL 

32vr 



t 2 j. 2 T ^0 

2^2 8M* 



z 02 



3/3 2 Ml 



( e p 2) (eV) 
A) i 2 



(51) 



where 7^ e p ^ and P ^ are pion-pion T-matrix elements in the isospin 
basis. The above three contributions have physically distinct origins; the 
first involves the direct effect of electromagnetism on the 1 = 2 decay ampli- 
tude, the second arises from final state scattering in which electromagnetism 
induces leakage from / = to / = 2, and the third is due to the shift in 
two-pion phase space produced by the electromagnetic mass shift pf . 



5 Final Results and Concluding Remarks 

Despite the presence of many unknown finite counterterms, it is possible 
to apply the numerical results of Eq. ( |50|) and obtain rough estimates of 
the EM corrections. The reasoning is that since the physical amplitudes 
are independent of the scale fj,, there must be compensating //-dependence 
between the chiral logarithms of Eq. (pi) and the counterterms. Therefore 
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the counterterms must be at least of the same order-of-magnitude as the 
chiral logs or even larger. We have adopted the operational procedure of 

(ct) 

assuming the counterterm contribution a\ vanishes at the scale fi = M p , 
and we assign an uncertainty given by ±|q|. This leads to the numerical 
values 

5(Al m ) = (0.024 ± 0.026) • 10" 7 M K o , 
5{Al m ) = (0.015 ± 0.022) • 10~ 7 M K o , 
5{A+ cm ) = (-0.005 ± 0.005) • 10~ 7 M R0 , (52) 
A 5/2 = (0.012 ± 0.016) • 10~ 7 M K0 , 

with A = (5.458 ± 0.012) -lO" 7 M K o and A 2 = (0.2454 ± 0.010) • 10" 7 M K o. 
Specifically, for the EM shift 6(A% em /A 2 ) calculated in Ref. |], we now have 
the extended result 

— '- = -(2.0 ±2.2) % . (53) 
A2 

If one allows for the uncertainty in g e mw in addition to those in the coun- 
terterm values, we find 



5(A 



2 ) 



(2.01^) % . (54) 



A 2 

In the numerical findings of Eqs. (|52|)-(|54[), the error bars are seen to be 
almost as large or larger than the signal. In our opinion, this is the best 
that one can do within a strict chiral perturbation theory approach. 
Our results illustrate several general features: 

1. Since the central values of the amplitudes have SA^ 171 7^ SA^ 171 ', the 
electromagnetic loop corrections are seen to produce AI = 5/2 effects, 
although the uncertainties of the counterterm values overwhelms the 
numerical result. 

2. A phenomenological analysis [|l?J] based on S- wave pion-pion scattering 
lengths and forward dispersion relations gives 5q — S 2 = (42 ±4)°. 
Yet an isospin analysis of K — > tttt decays yields 5o — 5 2 = (56.7 ± 
3.9)°. Presumably this difference of nearly 15° can be reconciled by 
subtracting EM effects from the K — > tttt decays. The main EM shift 
should be in 5 2 as only this angle experiences a AI = 1/2 enhancement. 
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Using Eq. (^T|) to calculate the angle 72 of Eq. (|3|) , we find 

Af p2) (3 



72 ,(eV) 32^ 



(e 2 p 0) 2^2 ( e p 2) 

02 3/3 2 M|. i2 



4.5° . (55) 



This evaluation, valid at order e 2 p°, is seen to worsen the discrepancy 
between the two determinations. To reveal the explanation behind 
this puzzle requires more work. || 

3. Finally, the most important implication of these estimates is that the 
electromagnetic shifts in A2 are not large, being only a few percent. 
Naive estimates allow the possibility that this shift could be much 
larger, perhaps even being a major portion of Ai- Our previous work 
at the leading order in the chiral expansion yielded a small effect. One 
motivation of the present calculation was to see if the next order effects 
upset this conclusion. Our estimates show that the natural size of the 
shift in A2 remains at the few percent level. 

This has been a complicated calculation with many different lagrangians, de- 
scribing different aspects of electromagnetic physics, required to obtain the 
full effect. These include explicit photon loops, mass shifts in the mesons 
propagating in loops and the short-distance electroweak interaction. The 
chiral power counting was crucial in sorting out which effects must be in- 
cluded for a consistent calculation. The resulting structure is universal and 
model independent. However, it is a prelude to more fully predictive ap- 
plications, as there remain unknown low energy constants which are not 
predicted by chiral symmetry alone. Different models can be used to es- 
timate the renormalized constants which appear in the chiral lagrangians, 
and these model predictions can then be readily translated into the physical 
amplitudes through the use of our calculation. In a following publication, 
we attempt to describe the extent that this may be accomplished using 
dispersive techniques to match long and short distance physics Jl. 
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